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ABSTRACT 

We report the detection of correlated anisotropics in the Cosmic Far-Infrared Background at 160 /im. 
We measure the power spectrum in the Spiteer/SWIRE Lockman Hole field. It reveals unambiguously 
a strong excess above cirrus and Poisson contributions, at spatial scales between 5 and 30 arcminutes, 
interpreted as the signature of infrared galaxy clustering. Using our model of infrared galaxy evolution 
we derive a linear bias h = 1.74 ± 0.16. It is a factor 2 higher than the bias measured for the local 
/i?^ 5" galaxies. Our model indicates that galaxies dominating the 160 /zm correlated anisotropics are 
at 2; ~ 1. This implies that infrared galaxies at high redshifts are biased tracers of mass, unlike in the 
local Universe. 

Subject headings: infrared: galaxies -galaxies: evolution - (cosmology:) large-scale structure of uni- 
verse 
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1. INTRODUCTION 

The discovery of the Cosmic Far-Infrared Background 
(CIB) in 1996, together with recent cosmological surveys 
from the mid-infrared to the millimeter has revolution- 
ized our view of star formation at high redshifts. It has 
become clear, in the last decade, that infrared galaxies 
contribute to a large part of the whole galaxy build- up 
in the Universe. Since the discovery of the CIB, new 
results on the identification of the sources contributing 
to the CIB, their redshift distribution, and their nature, 
are coming out at increasing speed, especially through 
mult i- wav elengt h analysis (see for a review Lagache et 
al. ()2005f )l. Stacking analysis are also very promising 
to p robe the CIB source popu lations (e.g. Dole et 
al. l[2006h . Wang et al. (I20M 1. Dye et al. l[200l V 
However, up to now, very little information is available 
on the clustering of infrared galaxies, although getting 
information on the clustering is essential to understand 
their formation process and to see how they relate to 
the other galaxy populations. 

The first three-dimensional quantitative measurements 
of the clustering strength of Ultra and Hyper Luminous 
Infrared Galaxies (ULIRGs, HyLIRGs) at hi gh red shifts 
(z > 1.5) have been made by Blain et al. (l2004f) . Far- 
rah et al. (|2006D and Magliocchetti et al. (|2007f) . These 
studies show that ULIRGs and HyLIRGs are associated 
with the most massive dark matter halos at high red- 
shifts, unlike in the local Universe where the star for- 
mation is quenched in the densest environments. CIB 
anisotropy observations provide a powerful complement 
to direct high angular resolution observations of individ- 
ual sources. CIB fluctuations measure, at large angu- 
lar scales, the linear clustering bias and, at small angu- 
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lar scales, the nonlinear cluste ring within a dark-matter 
halo (Cooray & Sheth ()2002D '). They thus probe both 
the dark-matter halo mass scale and the physics gov- 
erning the formation of infrared galaxies within a halo. 
However, up to now correlated anisotropics have never 
been firmly detected. In the far-infrared, detection of 
anisotropics is l imite d to the Poisson co ntribu tion (La- 
gache & Puget (I2000f). Matsuhara et al. (|2000D . Miville- 
Deschenes et al ()2002[ ))^. 

We report the detection of CIB correlated anisotropics at 
160 /im in the Spitzer SWIRE Lockman Hole field and 
give a first constraint on the linear bias. The paper is 
organised as follow: data are presented in Sect. [2l The 
power spectrum is analysed in Sect. [3l Finally a sum- 
mary and discussion are given in Sect. 3] 
Throughout this paper we use the cosmological parame- 
ters h 0.71, r^A = 0.73, ri„ = 0.27. For the dark mat- 
ter linear clustering we set the normalization to ag — 0.8. 

2. THE LOCKMAN HOLE SWIRE FIELD: MAP AND 
POWER SPECTRUM 

SWIRE has surveyed 49 square degrees distributed 
over 6 fields in th e northern and southern sky (Lons- 
dale et al. (|2004f )). The Lockman Hole is the largest 
field with the lowest cirrus emission. It covers about 10 
square degrees at 160 /im. 

2.1. Data reduction and map 

Raw data wer e redu ced using the Data Analysis Tool 
(Gordon et al. ^M)) version 2.71. We systematically 
removed each DCE (data collection event) after the stim- 
ulator flash to minimize the latency effect. We use the 
last calibration factor (44.7) to convert MIPS units to 
MJy/sr. Data were finally projected on a grid with 
15.95 arcsecond pixels. The map is shown on Fig. [TJ 
Further processings were necessary prior to measuring 
the power spectrum. We first needed to remove residual 
stripes. This has been done ex actly in the same way as 
in Miville-Deschene & Lagache ()2005f ). This paper shows 

^ Note however that Grossan &; Smoot 1)20061 ) report the detec- 
tion of the clustering signature at 160 /tm. 
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that our method efhciently removed residual stripes with- 
out affecting the astrophysical signal. We then remove 
all sources with Sieo >200 mJy (200 mJy is the high re- 
liability threshold, as detailed in Surace et al. (2005)). 
For this purpose we use DAOPHOT to detect the sources 
(the image was wavelet filtered prior to the detection). 
We then measure the fluxes using aperture photometry. 
After a fine centering on the sources, we integrate within 
25". Sky values are estimated in an [80"-110"] annu- 
lus. We compute the aperture correction - which is 2.02 
- using an effective instrumental function (that we call 
PSF for Point Spread Function) measured directly on a 
MIPS 160 /im map. Using an effective PSF rather than 
the PSF computed using the STinyTim program* is im- 
portant to take into account the survey strategy. In the 
Lockman Hole SWIRE field, the signal-to-noise ratio was 
not high enough to accuratly measure the PSF. We thus 
use the GTO/CDFS field in which the integration time 
is 6 times that in the Lockman hole SWIRE field. We 
checked that our measured fluxes at 160 /im were in very 
good agreement with the SWIRE DR2 catalog (better 
than 10% on average). On Fig [1] is shown the final map 
that will be used to compute the power spectrum. 

2.2. Power Spectra 

There are four contributors to the power spectrum at 
160 /im: cirrus emission, Poisson (shot) noise from dis- 
crete unresolved sources, CIB clustering (if any), and 
instrumental noise. If the noise and the signal are not 
correlated, the measured power spectrum P{k) follows: 

P{k)=-f{k)[P,,rrus{k)+Ps ources lus 

(fc)] +iV(fc) 

(1) 

where k is the 2D wavenumber (fc = \J k'^ + fc^, ex- 
pressed in arcmin ^), Pcirrus(k\ Psources^ Pclusik) 

and N{k) are respectively the power spectrum of the 
dust emission, the shot noise from unresolved sources 
(constant), the clustering and the noise. The factor 7(fc) 
represents the power spectrum of the PSF. To isolate 
the astrophysical components, we have to determine 
N{k) and7(fc). 

The noise power spectrum N{k) is computed by sub- 
tracting two maps of exactly the sa me region as detailed 
in Miville-Deschenes et al. ()2002f) . We construct two 
maps using the even and odd scans. As expected, N{k) 
and P{k) meet at small scales (fc ~ 1 arcmin"^) where 
the signal is noise-dominated. The noise power spectrum 
N(k) is subtracted from the raw power spectrum P{k). 

One of the critical issues is to correct the power spec- 
trum from the PSF 7(fc). The PSF at 160 /xm computed 
using the STinyTim program is very accurate but does 
not include any effect induced by the observing strategy. 
We have therefore also extracted direcly the PSF from 
the data (as discussed in Sect. [5]). The comparison 
of the power spectrum corrected by these two PSF is 
shown on Fig. [51 They are in very close agrement but 
we can notice that the effective PSF gives a better result 
(i.e. a flat power spectrum for 0.25 < k < 0.8 arcmin"^, 
as expected from Psources)- 

* http: / / ssc.spitzcr.caltcch.edu/archanaly/contributed/stinytim/ 



The error bars (shown in Fig. [3]) are estimated using 
a frequentist approach. Mock signal plus noise maps are 
generated and analysed with the same pipeline as for 
the data. This gives a set of power spectra of which we 
compute the covariance matrix. The diagonal elements 
of the covariance matrix are the errors on the measured 
power spectrum. 

3. POWER SPECTRUM ANALYSIS 

3.1. Adding low spatial frequency data to constrain the 
cirrus component 

Several studies show that the cirrus component dom- 
inate the power spectra at large scale s for k<0. 01 
arcmin-i (e.g. Miville-Deschenes et al. (|2002l : I2007D ). 
With the SWIRE data only, in such a low interstellar 
dust column density field, Pdrrus cannot be constrained. 
Larger maps are needed. We therefore compute the 
power spectrum of a large (--200 Sq. Deg.) IRIS/IRAS 
100 /im maps (Miville-Deschenes & Lagache (2005)). 
The SWIRE Lockman Hole field is embedded in this large 
IRIS map so that the average 100 /im dust emission is 
the same in the SWIRE and larger map (4% difference). 
Having the same average brightness is important since 
the normalisation of the cirrus power spectrum in the 
ver y diffu se region scales as I^qq (Miville-Deschenes et 
al. (|2007l )). We compute the power spectrum of the 100 
fim map after re movin g the bright sources as in Miville- 
Deschenes et al. ()2007[ ). We keep only the largest scales 
{k <9 10"'^ arcmin""'^), where we have only the contribu- 
tion from the cirrus component - the CIB being negligi- 
ble at these very large scales - and multiply the power 
spectrum by the average dust emission color (/i6o/-fioo)^- 
The color has been computed using DIRBE and FIRAS 
data. We compute the average |6| > 40° spectrum of the 
dust e mission correlated with the HI gas as in Lagache 
(|2003f) . We then fit the peak of the dust emission spec- 
trum to get the color. We obtain /i6o//ioo=2.06. If we 
take |6| > 30°, the color varies by ~10%. We show on 
Fig. [5] the 160 fim power spectrum derived from IRIS 
data together with the 160 /im MIPS power spectrum. 
The spectra agree impressively well. We can thus use 
this extended P{k) to constrain the cirrus contribution. 
In the following, the two spectra are stiched so that we 
use one spectrum from k^O.OOl to 1 arcmin^ ^. 

3.2. Detection of an excess at intermediate scales: 
signature of correlated anisotropics 

The cirrus component follows: 

f kY 

Pcirrus{k) = Pf) - I (2) 

\koJ 

where Pq is the power spectrum value at /cq^O.OI 
arcmin~^. Pq and (3 are determined by fitting 
the power spectrum (see Sect. 13. 3p . We obtain 
Pq = (2.98 ± 0.66) X lO^^ Jy^sr and /? = -2.89 ± 0.22. 
The normalisation can be converted at 100 fim 
using the 160/100 color given above; we obtain 
Po(100/im) — 7 X 10^ Jy^/sr. Considering that the mean 
cirrus value at 100 /im in our field is Iioo=0.51 MJy/sr, 
our measured PQ{lOOfim) and /? are in ex cellent agree- 
ment with Miville-Deschenes et al. (2007). The power 
spectrum of the cirrus component (Eq. ^ is displayed 
on Fig. [2] The measured power spectrum clearly has 
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an excess of power w.r.t. the cirrus contribution for 
k >0.3 arcmin"^. We interpret this strong excess as the 
signature of correlated CIB anisotropies. 

W e mod el the correlated anisotropies following Knox et 
al. ()2001[ )^. Using the three-dimensional, linear-theory 
power spectrum of dark matter density fluctuations to- 
day, Pm(^) the power spectrum of CIB anisotropies can 
be written as: 

= J ^^a'iz)-f{i^,z)b'PM{k)\k=i/rG'{z) ^ P,Uk) 

where r is the comoving proper-motion distance, k the 
3D wavenumber {k ~ ^/k'^ + k^ + k^, in Mpc~^), a(z) 



the scale factor, z) is the mean infrared galaxy emis- 
sivity per unit of comoving volume, and G{z) is the lin- 
ear theory growth function, i is the angular multipole 
sets using the Limber approximation, k —_t/r. We as- 
sume that the fluctuations in emissivity 5j / j are a biased 
tracer of those in the mass and introduce the bias param- 
eter h, that we assume independent of redshift and scale: 

v,z) <?p(k,z) 
^ = bx _ 4) 

where p is the dark matter density field. We compute 
the emissivity usi ng the infrared galaxy evolution model 
of Lagache et al. (l2004l ). This model, valid in the range 
3-1000 //m, is in very good agreement with mid-IR to 
far-IR number counts, CIB observations, resolved sources 
redshift distributions and local luminosity functio ns an d 
their evolution up to z ~ 2 (e.g. Lagache et al. ()2004D . 
Caputi et al. ((20061), Dole et al. (|2006f) . Frayer et al. 
dlOOe), Caputi et al. (|2007f) ). Fixing the cosmology, the 
only unkown parameter in Eq. [3] is the bias h. 

3.3. Measuring the bias 

We fit simultaneously Pq, P, b and Psources using the 
non-linear least-squares curve fitting mpfit program^. 
We obtain Pq = (2.98 ± 0.66) x 10^^ JyVsr, /3 = -2.89 ± 
0.22, b = 1.74 ± 0.16, and Psources = 9848 ± 120 JyVsr. 
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Pq and P have been discussed in Sect, 
agrees quite well with previous ISOPHOT determina- 
tion at ro ughly the same 5*160 threshold (Matsuhara et 
al. ()2000D V The fact that Pq, /3, and Psources are m very 
good agreement with previous measurements give us con- 
fidence in our measurement of the linear bias 6 1.7. It 
is well known that in the local Universe infrared galaxies 
are not bia sed tr acers of the mass. For example, Saun- 
ders et al. (|1992[ ) found &cr8=0.69±0.09 for IRAS galax- 
ies. Assuming cr8=0.8 gives b = 0.86. This bias is roughly 
comparable to the bias of the SDSS galaxies at z ~ 0.1 
(6 - 1.1, Tegmarket al (|2004D V We measure an average 
bias about 2 times higher in the CIB anisotropies. Figure 
3] shows the predicted redshift contribution to the cor- 
related anisotropies. At fc=0.05 arcmin"^, anisotropies 
from 0.7 < z < 1.5 infrared galaxies contribute for more 
'--^65%. Lower redshift galaxies contribute for less than 
5%. This shows that infrared galaxies at z~l are much 
more biased (~2 times) than locally. 



4. SUMMARY AND DISCUSSION 

We presented the power spectrum measured in the 
Spiteer/SWIRE Lockman Hole field at 160 /im. It is 
very well reproduced by the contribution from three 
components, cirrus, correlated CIB, and Poisson noise. 
The cirrus and Poisson contributions are very close to 
previous measurements. We measure the linear bias, 
b = 1.74 ± 0.16. This bias is likely to be that of infrared 
galaxies at z ~ 1 since z 1 galaxies are dominating 
the contribution to the correlated CIB anisotropies 
(Fig. HI and local galaxies are "anti-biased" . Such a 
bias is analog (but somewhat higher) to the bias of the 
red optical galaxy population at z ~ 1. For example, 
Marinoni et al. f|2005) measured a bias of 1.6 at z ~ 1.2 
for the red {B - I) > 1.5 galaxies in the VVDS. Blue 
galaxies at those redshifts are less biased with a relative 
bias between red and blue population of 1.4. 

The very strong evolution of the bias (from ~1.7 at 
z 1 to 0.86 at z=0) shows that as time progresses and 
the density field evolves, nonlinear peaks become less rare 
events and galaxy formation moves to lower-sigma peaks. 
Thus galaxies become less biased tracers of the mass den- 
sity field. Moreover, for the infrared galaxy population 
it is likely that an additional mecanism contributes sig- 
nificantly to the "debiasing" at low redshift. Galaxies 
in dense environments are found to have suppressed star 
formation rates (thus no or low infrared emission) and 
early morphological types compared with those in the 
field. Environmental effects in particular are important 
in quenching the star f ormati on through gas stripping 
(e.g. Postman et al. (j2005l )). though on some clus- 
ter outskirts, some star formation goes on (Due et al. 
([2002), Coia et al. (2005)). The high bias found for in- 
frared galaxies at z ~ 1 shows that star formation rates 
of galaxies are increasing wi th the environment (as also 
shown by e.g. Elbaz et al. l)2007f) ). The "merger bias" 
is an alternative, but somehow physically linked, way of 
boosting the bias at high redshift. Clustering of objects 
that have undergone recent mergers can be enhanced rel- 
ative to the clustering of individual halos o f comp arable 
masses (e. g. Fu rlanetto & Kamionkowski ()2005[ ): Wet- 
zel et al. (|2007l )). It is known that the star formation 
in infrared galaxies is triggered to some extent by merg- 
ers in dense environments (at z~l, 30-50% of luminous 
infrared galaxies are major mergers). Thus, the high 
measured bias may also point to "merger bias" . 
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^ Making a much more complex description of the correlated 



the scope of this paper 

clustering within the same dark matter halo) and the bias is beyond 



CIB anisotropics (as for example adding the contribution from the k , ; / v ■ ■ j ; ■ /-ji /■ n 1,4 i 

f \ f b u http://cow.physics.wisc.edu/~craigm/idl/idl.html 
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Fig. 1. — Left: 160 /^m SWIRE map after standard data reduction. Right: Final map used to compute the power spectrum. Residual 
stripes have been corrected as in Miville-Deschenes & Lagachc (.2005. ) . Sources with Sieo >200 mjy have been removed. 




Fig. 2. — Total power spectrum measured in the Lockman Hole at 160 /xm. In black (diamonds), the M/PS'/SWIRE 160 /xm power 
spectrum (P(k) — N(k), see Eq. 1), in purple (triangles) the IRIS/IRAS enlarged Lockman Hole field power spectrum at 100 fj,m, scaled 
to the power spectrum at 160 fim using the 160/100 dust color measured at high latitudes (|fe| >40°). In red and blue, the power spectra 

corrected from the STinyTim and measured PSF respectively. The green dashed-Sdotted line represents the best fit cirrus 

power spectrum, as computed in Sect. 13.31 
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Fig. 3. — Power spectrum in the Lockman Hole at 160 fim (black diamonds with error bars) with the three components: cirrus (green 
dashed-Sdotted line), clustering (blue, dashed-dotted line), Poisson (light blue, dashed line). The red continuous line is the sum of the 
three components. 
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Fig. 4. — Redshift contribution to the correlated anisotropics at 160 fim for b = 1.7. 



